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Introduction 

The  scientific  objective  of  this  research  is  to  explore  the  feasibility  of  assessing  tumor  status  during 
breast  can  cer  surgery  u  sing  t  argeted  f  luorescence-labeled  t  racers.  Under  normal  surgical  c  onditions, 
surgeons  rely  on  gross  examination  of  the  surgical  cavity  and  the  excised  tissue  to  help  guide  the  extent  of 
tissue  r emoval.  Post-surgical  analysis  us ing  thin  slices  of  tissue  stained  for  phy siochemical  features  is 
required  to  provide  definitive  assessments  of  the  distance  between  tumor  tissue  and  the  edge  of  the  excised 
tissue,  also  known  as  tumor  margin.  The  tumor  margin  status  has  a  significant  impact  on  the  risk  of  cancer 
recurrence,  and  a  diagnosis  of  “narrow”  or  “involved”  margins  requires  the  patient  to  return  for  additional 
surgical  resection.  Rates  of  repeat  surgeries  to  clear  tumor  tissue  are  surprisingly  high,  approaching  50% 
for  the  patient  population  undergoing  breast-conserving  surgery  for  ductal  carcinoma  in  situ  (DCIS),  a  type 
of  breast  cancer  that  has  yet  to  become  invasive.  We  hypothesize  that  targeted  fluorescent  drugs  coupled 
with  an  appropriate,  rapid  fluorescence  imaging  system  will  inform  the  surgeon  or  surgical  pathologist  of 
tumor  margin  status  in  excised  specimens  during  the  surgical  procedure.  Eliminating  the  need  to  inject  or 
otherwise  administer  the  drug  to  the  patient  would  allow  the  use  of  fluorophores  that  do  not  have  proven 
safety  profiles,  thus  removing  a  large  barrier  to  widespread  adoption  of  the  procedure  in  clinical  practice. 
The  project  involves  numerical  modeling  of  light  propagation  to  design  an  optimized  surgical  specimen 
scanner,  animal  experiments  to  investigate  the  behavior  of  different  targeted  fluorophores  applied  topically 
to  the  tissue  specimens,  and  a  pilot  study  with  human  tissue  in  the  later  stages  of  the  project.  Ultimately,  it 
is  our  hope  that  infonnation  provided  by  this  imaging  reduces  re-excision  rates. 

Body 

Progress  ha  s  be  en  m  ade  t  oward  t  he  pr  oposed  aims  of  this  project  including  de  velopment  of  light 
propagation  models  for  assessing  tumor  margin  status  (Task  3),  design  and  development  of  an  experimental 
imaging  system  (Task  3),  imaging  of  ex -vivo  tissues  from  animal  models  (Tasks  1  and  2),  and  modeling 
data  to  assess  sensitivity  and  specificity  metrics  of  imaging  in  vivo  (Task  4).  A  description  of  the  primary 
accomplishments  and  ongoing  efforts  follows. 

Light  propagation  modeling  and  diagnostic 

A  critical  aim  of  the  project  is  to  detennine  feasibility  of  detecting  tumor  margins  in  excised  tissue 
during  surgery.  Toward  that  end,  we  have  developed  sophisticated  light  modeling  algorithms  to  simulate 
the  light  pr  opagation  i  n  s  ections  of  t  issue  e  ncountered  dur  ing  s  urgery.  These  al  gorithms  ar  e  p  art  o  f  the 
NIRFAST  software  p  ackage  which  is  a  f  inite  e  lement-based  m  ethod  pa  ckage  ba  sed  on  t  he  di  ffusion 
approximation  of  photon  propagation  in  tissue.  Originally  developed  for  in  vivo  imaging  in  relatively  large 
tissue  volumes,  this  software  has  been  modified  to  investigate  reflectance-mode  fluorescence  imaging  for 
surgical  settings  as  part  of  this  project.  Accurate  light  modeling  of  tissue  sections  allows  optimization  of 
the  intra-surgical  imaging  system. 

Figure  1  illustrates  the  modeling  capabilities  developed.  A  simulated  test  domain  was  creating  from 
an  a  rchive  phot  ograph  of  a  br  east  t  issue  s  ection  (Fig.  1(a)).  By  e  stimating  a  tis  sue  th  ickness,  the  2-D 
outline  of  the  tissue  from  the  photograph  was  used  to  generate  a  3-D  finite  element  mesh  of  46,000  nodes. 
Four  simulated  tumor  fragments  were  added  to  the  simulated  tissue  volume  at  different  depths  from  the 
tissue  s  urface.  S  emi-transparent  surface  renderings  containing  i  nterior  c  ross  s  ections  i  n  t  wo  di  fferent 
planes  to  show  the  location  and  shapes  of  the  simulated  tumor  regions  are  presented  in  Figs.  1(b)  and  (c). 
In  this  example,  tumor-to-background  fluorescence  intensity  was  assumed  to  be  20: 1  for  all  four  fragments. 
Figure  2  (d)  shows  the  s  imulated  illu  mination  in  tensity  o  f  th  e  e  xcitation  s  ource  o  n  th  e  tis  sue  s  urface, 
modeled  as  a  dense  grid  of  point  sources.  This  excitation  illumination  induces  fluorescence  throughout  the 
tissue  volume.  The  remitted  fluorescence  is  shown  in  Fig.  2(e),  and  represents  what  the  surgeon  would  see 
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under  these  conditions.  M  ost  tumor  fragments  generate  substantial  fluorescence  contrast  on  the  surface. 
The  small  fragment  buried  further  in  the  tissue  is  not  detectable. 


Figure  1 .  An  archive  photo  of  a  t  issue  s  pecimen  i  n  (  a)  w  as  u  sed  tog  enerate  a  f  inite  el  ement  m  esh 
containing  simulated  tumor  fragments  shown  in  (b)  and  (c).  T  he  remitted  fluorescence  from  the  tumor 
tissue  excited  by  an  illumination  source  shown  in  (d)  is  presented  in  (e). 

Spectrally  resolved  wide-field  imaging  for  assessing  tumor  depth 

The  ability  to  quickly  determine  tumor  depth  from  the  tissue  surface  either  in  a  fresh  specimen  or  in 
a  surgical  cavity  may  benefit  efforts  to  reduce  patient  callback  after  breast  surgery.  One  approach  towards 
quantify  tumor  depth  is  to  take  advantage  of  the  spectral  distortion  of  the  fluorescence  peak  as  it  propagates 
through  tissue  (1-3).  This  is  due  to  the  spectral  variation  of  optical  properties.  Fluorescence  from  deeper 
tumors  will  be  distorted  to  a  greater  degree  than  those  near  the  surface,  and  it  may  be  possible  to  exploit 
this  phenomenon  to  provide  an  estimate  of  depth  that  is  clinically  relevant. 

The  simulated  experiment  in  Fig.  2  demonstrates  the  spectral  distortion  of  an  embedded  tumor.  The 
large  gray  r  egion  r  epresents  a  2  -D  c  ross  s  ection  of  t  issue  w  ith  t  he  s  ources  a  ligned  a  long  the  t  op  of  t  he 
region.  These  sources  were  assumed  to  be  illuminated  simultaneously  to  represent  a  broad-beam  imaging 
system.  A  simulated  tumor  with  10: 1  contrastto  the  background  is  shown  at  a  relatively  deep  position 
within  t  he  t  issue.  I  n  t  his  s  tudy,  t  he  de  pth  of  t  his  t  umor  was  varied  over  a  r  ange  of  a  bout  1 8m  m,  as 
illustrated  in  the  figure.  At  each  depth,  the  spectrum  of  the  emitted  fluorescence  was  calculated  for  discrete 
positions  on  t  he  tissue  surface.  T  hese  spectra  are  shown  in  the  figure.  Obviously,  spectra  far  from  the 
position  of  t  he  t  umor  a  re  not  affected  n  oticeably  b  y  t  he  eh  anging  d  epth.  N  ear  t  he  t  umor  the  o  verall 
intensity  of  the  remitted  fluorescence  changes  dramatically  with  depth  and  there  is  a  noticeable  change  in 
the  shape  of  the  spectrum. 
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O  Source 
x  Detector 


Figure  2.  Three  fluorescence  spectra  remitted  from  the  tissue  at  different  spatial  points  are  shown  above 
the  blue  arrows.  The  depth  of  the  simulated  tumor  mass  was  varied,  as  can  be  seen  in  the  cross-sectional 
view  of  the  tissue  in  gray.  For  points  on  the  tissue  surface  near  the  tumor,  the  remitted  fluorescence  spectra 
change  noticeably  with  changing  tumor  depth 

A  simple  ratio  between  two  wavelengths  may  provide  depth  information.  Figure  3  shows  the  same  domain, 
again  illustrating  different  positions  of  the  same  tumor  as  in  Fig.  2,  above  which  is  reported  the  ratio  of  two 
wavelengths  in  the  emission  spectrum  resolved  spatially.  C  hanges  in  the  ratio  correspond  to  changes  in 
depth,  how  ever,  t  he  r  atiometric  m  easurement  i  s  not  uni  que  w  ith  d  epth.  O  ther  f  actors,  s  uch  as  r  elative 
position  on  t  he  t  issue  s  urface,  s  ize  of  t  he  t  umor,  a  nd  c  ontrast  w  ill  a  11  pi  ay  a  r  ole  i  n  t  he  r  atio  va  lue. 
However,  if  a  few  of  these  parameters  can  be  assumed,  which  is  feasible,  it  may  be  possible  to  decrease  the 
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non-uniqueness  of  the  problem.  This  is  part  of  ongoing  work  on  this  project. 


Tissue 

surface 


Figure  3.  The  ratios  of  two  wavelengths  in  the  e  mission  s  pectrum  of  the  fluorophore  a  re  plotted  for 
different  points  on  the  tissue  surface.  The  white  blotches  within  the  black  diagram  illustrate  the  different 
tumor  positions.  Detector  positions  on  the  tissue  surface  correspond  to  the  ratio  plot  directly  above. 


We  have  also  developed  an  innovative  mathematical  technique  to  estimate  depth  based  on  analytical 
solutions  t  o  t  he  di  ffusion  e  quation.  T  his  simplification  i  s  a  Iso  ba  sed  on  ratiometric  m  easurements  at 
different  wavelengths  but  has  significant  computational  efficiency  advantages  over  numerical  approaches 
and  therefore  can  be  integrated  easily  and  calculated  at  higher  rates  than  video.  Please  refer  to  the  attached 
paper  for  details.  While  the  paper  focuses  on  a  specific  organ  region  (brain  tumor  surgery  in  this  example), 
the  technique  is  generic  in  its  application  and  can  readily  be  applied  to  any  organ. 

Finally,  we  are  currently  working  to  integrate  more  accurate  light  propagation  models,  such  as  Monte 
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Carlo  simulations  for  irregular  domains,  into  the  software  package.  T  hese  models  will  provide  a  more 
accurate  representation  of  what  the  surgeon  will  see  with  an  intra-surgical  imaging  system  and  thus  allow 
further  refinement  of  these  system  configurations. 


Imaging  system  design 

An  initial  experimental  assessment  of  the  analytical  ratiometric  method  for  determining  tumor  depth 
has  been  performed  by  designing  an  instrument  for  wide-field  acquisition  of  multi-spectral  measurements, 
as  s  hown  i  n  F  ig.  4 .  In  this  e  xperiment,  fluorescence  ex  citation  was  done  us  ing  an  expanded  b  earn  o  f  a 
diode  1  aser  w  ith  w  avelength  Aex  =  635  nm  (Power  T  echnology  I  nc.,  L  ittle  R  ock,  A  rkansas).  R  emitted 
fluorescence  was  measured  at  five  wavelengths,  X  =  (675,  685,  695,  705,  715)  nm,  using  10  nm-wide  band¬ 
pass  interference  filters  (Omega  Optical,  Brattleboro,  Vermont)  inserted  into  a  motorized  wheel  coupled  to 
an  EXi  Aqua  charge-coupled  device  (CCD)  (Qlmaging,  Surrey,  British  Columbia).  The  camera  was  cooled 
at  0°C  and  has  a  full  well-depth  of  18,000  e-,  a  14  bits  digitizer,  a  dark  current  of  0.15  e-/pixel/seconds  and 
a  read-noise  of  6.5  e-  when  operated  at  10MHz  readout  frequency. 


CCD 

camera 


Broad  beam 
source 


(a) 


Excitation  light 


uuuuuu 


Figure  4.  (a)  N  on-contact  mu  lti-spectral  imaging  instrument  us  ed  f  or  f  luorescence-ratio  de  tection.  ( b) 

Light  transport  is  modeled  assuming  the  cavity  can  be  approximated  as  a  semi-infinite  diffusive  medium 
with  homogeneous  optical  properties. 


This  system  was  used  to  test  the  analytical  ratiometric  methodology  by  imaging  a  tissue-simulating 
blood  phantom  with  protoporphyrin  IX  as  a  fluorophore.  A  container  (dimensions:  150  mm  x  75  mm  x  61 
mm)  was  filled  with  water  containing  a  v  olume  f  raction  o  f  1  %  in  tralipid  (  Baxter  H  ealthcare  C  orp., 
Deerfield,  Illinois)  and  2%  porcine  blood.  An  optically  transparent  cylindrical  inclusion  (diameter:  5.5  mm, 
height:  10  mm)  was  filled  with  20  jug/ml  of  PpIX  (Frontier  Scientific,  Logan,  Utah)  dissolved  in  dimethyl 
sulfoxide  (DMSO).  A  moving  stage  allowed  the  simulated  tumor  to  be  immersed  at  different  depths.  The 
medium  was  illuminated  using  an  expanded  laser  beam  (diameter:  13  m  m,  P  =  5.5  m  W,  /.cx  =  635  nm  ). 
Fluorescence  images  at  five  different  wavelengths  were  acquired  (d  =  1  to  15  mm).  A  PpIX  spectrum  was 
acquired  to  be  used  as  a  reference  that  was  not  corrupted  by  light  propagation  in  the  medium.  Associated 
with  each  fluorescence  image,  a  p  re-contrast  image  (inclusion  filled  with  DMSO  only)  was  acquired  and 
subtracted  from  the  c  orresponding  i  mage  t  o  a  ccount  f or  e  xcitation  1  ight  filter  bl  eed- through.  Fig.  5(b) 


shows  the  measured  ratios  for  two  sets  of  wavelengths:  (A,i,  A2)  =  (675  nm,  685  nm)  (circles)  and  (/.|,  A,2)  = 
(675  nm,  695  nm)  (squares).  The  dotted  lines  correspond  to  simulated  ratios  computed  with  the  analytical 
approximation  using  the  optical  properties  experimentally  retrieved  with  a  fiber-optics  probe  (Fig.  5(a)). 
Qualitative  inspection  of  the  data  suggests  a  reasonable  match  between  theory  and  experiment,  although  a 
significant  level  of  noise  is  observed  mainly  due  to  excitation  light  filter  bleed-through,  limited  light  power 
delivered  to  the  tissue  surface  and  variations  of  t  he  opt  ical  pr  operties  of  t  he  bl  ood  pha  ntom  dur  ing  t  he 
experiment  (as  much  as  20%  for  the  absorption  coefficient).  Using  a  ( non-parametric)  Speannan’s  rank 
correlation  te  st  a  s  ignificant  co  rrelation  w  as  f  ound  b  etween  t  he  ex  perimental  d  ata  an  d  t  he  t  heoretical 
predictions.  The  correlation  coefficients  and  P-values  obtained  were  rs  =  0.90,  P  =  0.0046  (two-tailed)  for 
(A-i,  A,2)  =  (675  nm,  685  nm)  and  rs  =  0.83,  P  =  0.015  (two-tailed)  for  (A-i,  A,?)  =  (675  nm,  695  nm).  A  linear 
regression  was  performed  on  the  data  with  goodness-of-fit  parameters  r2  =  0.84  and  r2  =  0.64  for  (A,i,  7,?)  = 
(675  nm,  685  nm)  and  (7,i,  7,2)  =  (675  nm,  695  nm),  respectively.  Another  statistical  analysis  (F-test)  was 
performed  to  test  for  the  presence  of  a  non-zero  slope  in  the  relationship  between  T  and  d.  We  found  that 
for  both  sets  of  wavelengths  the  data  was  significantly  different  from  a  horizontal  line  (P<0.05). 


T - 1 - 1- 

5  10  15 

depth  (mm) 


—  705  -  675  nm  (theory) 

—  685  -  675  nm  (theory) 
•  705  -  675  nm  (meas) 
■  685  -  675  nm  (meas) 


(a)  (b) 

Figure  5.  Experimental  data  acquired  using  a  tissue-simulating  phantom:  (a)  Optical  properties  (absorption 
and  r  educed  s  cattering)  m  easured  w  ith  t  he  a  fiber  probe,  ( b)  S  imulated  ( theory)  v  ersus  m  easured  (meas) 
fluorescence  ratios  for  two  sets  of  wavelengths  showing  correlation  with  depth. 


Animal  experiments 

In  c  ollaboration  w  ith  D  r.  Hoopes’s  lab,  excised  tissue  from  a  small  number  (N=4)  of  animals  was 
imaged  on  the  Licro  Odyssey  specimen  scanner.  A  fter  tumors  reached  approximately  7  mm  across,  the 
animals  were  sacrificed  and  the  tumor  immediately  removed  with  adjacent  nonnal  tissue.  All  tissue  were 
soaked  for  5-20  minutes  in  a  solution  of  Licor  IRDye800CW-2DG  and  then  rinsed  for  5-20  minutes  in  a 
PBS  s  olution.  T  he  t  issues  w  ere  t  hen  s  canned  on  t  he  fluorescence  s  canner  a  nd  w  hite  1  ight  phot  os  t  aken 
without  changing  the  geometry  of  the  specimens.  F  igure  6  shows  the  white  light  and  fluorescence  image 
pairs  for  specimens  taken  from  two  animals.  T  he  two  central  images  in  (d)  are  of  a  cross-section  of  the 
tumor  slices  in  two  after  being  soaked  and  rinsed.  These  images  indicate  that  the  drug  does  not  appear  to 
penetrate  the  tumor  well  and  normal  tissue  uptake,  or  auto  fluorescence,  seems  to  be  significant  with  this 
drug.  These  preliminary  results  suggest  that  this  optical  probe  may  not  have  the  specificity  required  for 
margin  detection;  however,  further  examination  is  required  to  be  conclusive. 
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Figure  6.  White  light  images  and  corresponding  fluorescence  images  for  tumors  excised  from  two  animals 
((a)  and  (b)  for  one  animal  and  (c)  and  (d)  for  the  second).  In  (c),  the  central  two  specimens  are  the  same 
tumor  sliced  down  the  middle  after  soaking  in  Licor  IRDye800CW-2DG,  showing  a  cross  section  of  the 
tumor  mass.  These  images  show  a  significant  amount  of  non-specific  uptake  of  the  probe  in  normal  tissue 
and  poor  penetration  into  the  tumor  mass. 

Working  with  Dr.  Hoopes,  Dr.  Julia  O  ’hara  and  their  colleagues,  I  have  gained  experience  in  the 
methods  used  to  implant  breast  cancer  tumor  lines  in  the  mammary  fat  pad  of  nude  mice.  As  the  study  has 
progressed,  we  have  switched  from  the  S  KBR-3  breast  cancer  line  to  a  HER2/neu  positive  B  T474  and 
MDA23 1  (HER2/neu  negative)  which  we  have  found  grow  more  readily  in  vivo.  The  use  of  a  n  egative 
control  is  an  important  addition  to  the  study  design.  Currently,  we  have  40  animal  subjects  growing  breast 
cancer  t umor  m odels  for  this  project,  20  w  ith  HER2/neu  positive  B  T474  tumors  and  20withMDA231 
cells.  Implantations  were  perfonned  with  5  x  106  cells  in  matrigel.  Estrogen  pellets  were  also  implanted 
under  the  skin  of  mice  inoculated  with  BT474  cells  to  promote  tumor  growth.  Tumors  in  these  populations 
should  be  mature  within  1-2  months  from  the  date  of  this  report. 

Assessment  of  diagnostic  performance  of  optical  imaging  of  breast  cancer  in  vivo 

One  component  of  the  project  involves  in  vivo  assessment  of  breast  using  MRI-guided  fluorescence 
molecular  tomography.  In  this  paradigm,  tissue  region  segmented  out  of  the  simultaneously  acquired  MR 
images  is  used  to  g  uide  t  he  d  iffuse  o  ptical  i  mage  r  econstruction  pr  ocess.  While  M  Rl-coupled  opt  ical 
tomography  has  matured  over  the  past  decade,  the  benefits  of  incorporating  MRI  structural  information  in 
the  image  recovery  of  optical  properties  has  not  been  quantified  beyond  anecdotal  case  studies.  In  this 
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study,  t  he  i  mprovements  of  M  RI  pr  iors  ove  r  s  tandard  di  ffuse  i  maging  without  spatial  guidance  are 
quantified  in  terms  of  diagnostic  performance.  A  simulated  breast  volume  made  up  of  fatty,  fibro-glandular 
and  tumor  regions  was  used  as  the  test  domain.  A  range  of  tumor  sizes  and  contrasts,  including  MRI  false 
positive  cases,  were  simulated  and  images  were  recovered  with  or  without  the  use  of  the  MRI  spatial  prior 
information.  Receiver  Operating  C  haracteristic  (  ROC)  an  alysis  w  as  u  sed  t  o  q  uantify  d  iagnostic 
performance  of  each  technique.  The  results  are  shown  in  Fig.  7  for  ROC  analysis  perfonned  without  MRI 
guidance,  with  MRI  priors  assuming  that  the  internal  structure  is  accurate,  and  with  MRI  priors  assuming 
that  the  internal  structure  of  the  fibro-glandular  tissue  is  incorrectly  dilated  by  1  pixel,  for  (a),  (b),  and  (c), 
respectively.  Area-under-the-curve  (AUC)  values  are  0.73,  0.85,  and  0.75,  respectively,  indicating  that  the 
use  of  MRI  priors  i  mproves  di  agnostic  pe  rformance,  pr  ovided  t  he  i  ntemal  s  tructure  obtained  f  rom  t  he 
modality  is  accurate.  Further  analysis  is  ongoing. 


Figure  7.  ROC  c  urves  f  or  di  agnostic  pe  rformance  of  s  tandard  di  ffuse  t  omography,  (  a),  i  mage-guided 
diffuse  t  omography  a  ssuming  a  ccurate  s  tructural  pr  iors,  (b),  and  image-guided  di  ffuse  t  omography 
assuming  inaccurate  structural  priors,  (c). 


Key  research  accomplishments 

•  Developed  m  odeling  a  lgorithms  for  s  imulating  ima  ging  o  f  fresh,  excised  t  issue  and  t  he  s  urgical 
cavity. 

•  Developed  a  method  for  identifying  depth  of  tumors  from  the  tissue  surface  in  a  s  urgical  setting. 
This  technique  is  applicable  to  tissue  volumes  in  the  surgical  cavity  or  fresh,  excised  tissue. 

•  Contributed  to  the  design,  construction,  and  experimental  validation  of  a  s  pectrally-resolved  intra- 
surgical  imaging  system. 

•  Acquired  pr  eliminary  da  ta  in  a  s  mall  n  umber  o  f  a  nimals  to  demonstrate  penetration  of  targeted 
fluorescence  p  robe  i  n  f  resh  ex  cised  t  issue.  Initial  r  eadings  do  not  s  how  a  hi  gh  1  evel  of  specificity  or 
penetration  depth  with  this  optical  probe. 

•  Assessed  diagnostic  performance  of  MRI-guided  in  vivo  optical  imaging  using  Receiver  Operating 
Characteristic  analysis. 
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Reportable  outcomes 

The  first  year  of  this  traineeship  grant  has  led  to  one  second  author  publication  manuscript  currently 
under  review  (Journal  of  Biomedical  Optics),  one  second  author  paper  in  preparation,  and  one  first  author 
paper  i  n  pr  eparation.  In  a  ddition  t  o  a  ttending  s  eminars  a  nd  c  onferences,  I  also  served  a  s  a  w  orkshop 
instructor  for  the  Workshop  on  Diffuse  Optical  Tomography:  NIRFAST  software  using  MATLAB  at  OSA 
Biomedical  O  ptics  ( April  12-14,2010),  and  an  organizer  and  workshop  i  nstructor  a  1 1  he  Image-Guided 
Spectroscopy  Symposium  and  Workshop  at  the  Thayer  School  of  Engineering,  Dartmouth  College  (July  19 
-23,2010). 

Conclusions 

Substantial  progress  has  been  made  on  the  numerical  modeling  and  system  development  aspects  of 
the  pr  oject.  T  hree  di  mensional  m  odeling  of  t  umors  i  n  e  xcised  t  issue  ha  s  be  en  a  ccomplished  a  nd  t  he 
introduction  of  spectral  resolution  in  the  imaging  strategy  shows  promise  for  quantifying  tumor  depth  in  the 
tissue.  A  pr  ototype  i  maging  s  ystem  ha  s  be  en  de  veloped  for  intra-surgical  ma  rgin  d  etection  w  ith  in  itial 
studies  de  monstrating  t he  potential  t o  de termine  t he  de pth  of  e mbedded  tumor  f ragments  i  n  t  he  e xcise 
tissue  or  surgical  cavity.  Pilot  animal  data  with  the  Licor  IRDye800CW-2DG  optical  probe  in  MDA23 1 
tumors  showed  poor  s  pecificity  a  nd  t  umor  pe  netration.  T  hese  results  are  not  promising  for  diagnosing 
tumor  margins  with  this  fluorescent  compound;  however,  further  investigation  with  this  probe  is  ongoing. 
While  the  animal  experiments  have  progressed  at  a  s  lower  pace  than  originally  proposed,  I  have  received 
substantial  training  in  this  area  and  with  40  animals  currently  growing  xenograft  tumors  in  the  mammary 
fat  pad,  the  pace  of  experiments  will  increase  sharply  within  the  next  1-2  months.  These  studies  will  also 
involve  the  trastuzumab-bound  NIR  fluorophore  which  is  expected  to  provide  higher  specificity  between 
tumor  and  normal  tissue  in  the  HER-2  positive  tumor  line.  Finally,  ROC  analysis  was  used  to  quantify 
imaging  performance  of  an  in  vivo  MRI-coupled  diffuse  tomography  imaging  system  and  these  data  are 
being  prepared  for  publication. 
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The  relevant  journal  article  manuscript  has  been  attached.  This  manuscript  is  under  review. 
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Abstract:  Analytical  solutions  to  diffuse  light  transport  in  biological  tissue  are  derived 
providing  a  closed  form  expression  to  estimate  the  depth  of  fluorescent  molecules 
based  on  the  spectral  deformation  of  diffused  near-infrared  measurements.  It  is  shown 
that  the  predicted  depth  is  linearly  related  to  the  logarithm  of  the  ratio  of  fluorescence  at 
two  different  wavelengths,  as  long  as  the  depth  is  beyond  a  few  millimeters  and  the 
tissue  is  relatively  homogenous.  In  agreement  with  theoretical  predictions,  experimental 
data,  acquired  on  a  tissue-simulating  blood  phantom  with  a  broad-beam  non-contact 
multi-spectral  imaging  system,  is  used  to  confirm  that  the  ratio  of  fluorescence  at  two 
different  wavelengths  correlates  with  depth.  Further  elucidation  of  the  relationship 
between  depth  and  deformation  of  fluorescence  spectra  is  provided  by  demonstrating 
that  that  the  slope  and  intercept  values  of  the  relation  between  fluorescence  ratio  and 
depth  is  a  function  of  tissue  absorption  and  reduced  scattering.  Those  findings  are  used 
as  a  basis  for  proposing  that  fluorescence  contrast  depth  assessment  can  be  achieved 
in  fluorescence-guided  neurosurgery  allowing  improved  intra-operative  delineation  of 
tumor  margins.  In  principle,  this  approach  could  also  be  applied  in  other  preclinical  or 
clinical  imaging  where  luminescence  contrast  depth  assessment  is  required  in  a  multi¬ 
pixel  sub-surface  imaging  geometry. 

Keywords:  fluorescence,  multi-spectral  imaging,  diffusion  theory,  image-guided 
surgery,  protoporphyrin  IX,  brain  tumors 
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1  Introduction 


Over  the  course  of  the  last  decade,  there  has  been  increased  interest  in  the 
development  of  pre-clinical  and  clinical  methods  allowing  tissue  resection  based  on 
intra-operative  fluorescence  imaging.  The  motivation  for  this  type  of  research  is  the 
application  of  optical  methods  in  the  operating  room  that  provide  surgeons  with 
pertinent  tissue  information  that  is  complementary  to  that  conveyed  by  conventional 
white-light  images.  A  clinical  application  of  fluorescence-guided  surgery  that  is 
becoming  more  widespread  consists  in  the  resection  of  brain  tumors  using  the  intrinsic 
metabolic  and  structural  changes  induced  by  the  administration  of  S-aminolevulenic  acid 
(ALA).  In  this  case,  exogenous  administration  of  ALA  overloads  the  heme  cellular 
pathway  leading  to  selective  accumulation  of  the  fluorescent  heme  precursor, 
protoporphyrin  IX  (PplX),  in  neoplastic  tissues.  Accumulated  levels  of  PplX  in  abnormal 
tissues  are  sufficiently  high  that  certain  types  of  tumors  can  be  visualized  through  a 
surgical  microscope  modified  with  a  fluorescence  imaging  channel  [1, 2], 

The  current  state-of-the-art  methodology  for  fluorescence-guided  brain  tumor 
resection  consists  in  using  a  modified  surgical  microscope  to  excite  the  PplX  molecules 
using  a  broad-beam  blue  light  source.  Photodetection  of  the  re-emitted  fluorescence  is 
then  achieved  using  a  charged-coupled  device  (CCD)  camera  coupled  to  a  high-pass 
filter  [3],  However,  an  inherent  limitation  of  this  in  vivo  fluorescence  imaging  method  is 
the  strong  light  absorption  in  the  visible  part  of  the  electromagnetic  spectrum,  typically 
making  measurements  sensitive  only  to  biological  events  occurring  within  a  few  layers 
of  cells.  Therefore,  the  current  state-of-the-art  for  fluorescence-guided  resection  is 
severely  limited  in  its  ability  to  identify  sub-surface  tumor  remnants.  Light  in  the  NIR, 
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however,  can  travel  much  deeper  into  tissue  because  there  tissue  absorption  decreases 
by  two  to  three  orders  of  magnitude,  when  compared  to  visible  light.  In  the  NIR  range, 
scattering  becomes  the  dominant  light  attenuation  mechanism.  Despite  the  resulting 
increased  depth  sensitivity  however,  planar  epi-illumination  NIR  images  are  difficult  to 
interpret  because  the  information  they  convey  is  surface-weighted  with  a  strong  bias 
toward  fluorescence  near  the  surface  [4], 

Moving  beyond  simple  planar  images  by  providing  surgeons  with  complementary 
information  relating  to  the  depth  of  fluorescent  sources  is  an  important  aspect  of  future 
developments  in  fluorescence-guided  surgery.  A  possibility,  when  it  comes  to  retrieving 
depth  information,  consists  in  developing  an  intra-operative  diffuse  fluorescence 
tomography  instrument.  Such  a  method  would  involve  a  laser  beam  scanned  across  the 
region-of-interest  (ROI)  on  the  tissue  surface.  For  each  illumination  point,  a  CCD  image 
would  be  acquired  providing  a  high  density  of  measurements  tomography  dataset  [5,  6], 
Given  proper  light  transport  modeling  can  be  achieved,  an  inverse  problem  would  then 
be  solved  providing  the  surgeon  with  intra-operative  tomography  images  with  tissue 
information  up  to  1-2  cm  underneath  the  surface  of  the  surgical  cavity.  However,  there 
are  several  impediments  to  the  integration  of  such  an  approach  in  neurosurgery.  In 
order  to  be  practical,  a  depth-resolved  approach  should  be  implemented  in  such  a  way 
that  the  surgical  workflow  is  minimally  disturbed.  From  that  point  of  view,  intra-operative 
diffuse  fluorescence  tomography  might  not  be  an  optimal  choice  because  of  the  large 
time  required  to  acquire  a  full  dataset  with  acceptable  signal-to-noise  characteristics,  as 
well  as  the  potentially  large  time  required  to  actually  reconstruct  and  properly  analyze 
the  tomography  images. 
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This  paper  presents  preliminary  results  in  the  development  of  a  depth-resolved 
fluorescence  imaging  approach  compatible  with  the  neurosurgical  workflow.  The 
proposed  methodology  is  a  novel  non-tomographic  approach  potentially  allowing  pixel- 
by-pixel  fluorophore  depth  estimation  to  be  obtained  at  an  imaging  rate  that  is 
sufficiently  close  to  real-time  and  that  introduces  minimal  computational  burden.  The 
proposed  methodology  consists  in  acquiring  two  fluorescence  images  at  different 
wavelengths  and  using  their  ratio  to  estimate  the  depth  associated  with  the  sources  of 
fluorescence  located  underneath  the  tissue  surface.  Studies  have  shown  that 
computational  depth  estimation  based  on  multi-spectral  data  is  feasible  [7,  8],  However, 
the  results  derived  from  those  studies  are  mainly  applicable  to  contact  point 
fluorescence  measurements  using  fiber-optics,  an  approach  that  is  not  compatible  with 
routine  use  in  the  scope  of  fluorescence-guided  surgical  applications.  In  comparison, 
the  developments  presented  here  provide  evidence  that  analytical  solutions  to  diffuse 
light  transport  can  be  used  to  estimate  depth  using  a  non-contact  imaging  system. 
Moreover,  the  depth-retrieval  formalism  developed  here  is  argued  to  be  amenable  to 
wide-field  CCD-based  photo-detection  in  a  manner  that  is  compatible  with  the  surgical 
workflow.  On  a  more  conceptual  level,  the  work  presented  here  fills  a  gap  in  the 
literature  by  providing  closed  form  expressions  from  which  the  relation  between  spectral 
deformations,  depth  of  fluorophores  and  intrinsic  optical  properties  can  be  formally 
understood. 

2  Methods  and  Results 

2.1  The  relationship  between  fluorescence  ratio  and  depth 
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An  example  of  a  non-contact  configuration  for  intra-operative  multi-spectral  imaging  is 
illustrated  in  Fig.  1(a)  where  the  optical  light  paths  can  be  made  compatible  with  those 
of  existing  surgical  microscopes.  In  this  configuration,  the  fluorescence  excitation  is 
generated  using  a  broad-beam  source  in  the  NIR  at  the  excitation  wavelength,  Xex. 
Remitted  fluorescence  signals  are  measured  at  two  wavelengths,  A,i,2>Xex,  in  the 
fluorescence  emission  spectrum  using  band-pass  filters  and  a  cooled  CCD  chip. 
Assuming  that  tissue  absorption  due  to  the  fluorophores  is  much  smaller  than  that  of 
other  native  chromophores  ( e.g .,  oxy-hemoglobin,  deoxy-hemoglobin)  the  light  signals 
at  each  emission  wavelength  can  be  modeled  using  the  expression 

Yem (R, X* , X)  *  QFe? (X)  J d3r  y/** (i r ) CF (?) Gem (R,r,SA,DA),  ( 1 ) 

Q 


where  QF  is  the  quantum  yield,  sFem(^)  the  emission  spectrum  and  CF(r)  the 
concentration  of  fluorophores  at  location  r.  Fig.  1(b)  shows  the  main  geometrical 
variables  used  in  Eq.  (1),  in  which  light  transport  is  modeled  as  a  diffusive  process.  The 
fluence  function  y/ex  is  the  excitation  light  field  and  Gem  is  the  diffusion  equation  Green’s 
function,  which  corresponds  to  the  radiant  exposure  in  response  to  a  light  impulse  at  r. 
For  boundary  conditions  associated  with  an  infinite  homogenous  medium,  it  is  given  by 


G™(R,r  ,8A  ,DX)  = 


exp(-|i?-r  1 1 8a) 
AkDa\R-?\ 


(2) 
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where  the  diffusion  constant  is  Dx  =  l/3(//^  +  ju*')  and  the  penetration  depth  is 

8‘  =  [9]- 

A  closed  form  expression  is  now  derived  from  which  a  depth  value  can  be 
attributed  to  each  point  imaged  on  the  tissue  surface.  This  is  performed  making  the 
approximation  that  all  of  the  fluorescence  emanating  from  the  surface  is  from  a  point¬ 
like  distribution  with  molar  concentration  CF  at  rs.  This  is  equivalent  to  setting 
CF(r)  =  CF8°\r  -rs)  in  Eq.  (1),  which  then  takes  the  form 


¥e,\R,r,X)^CFQFsTW^(K)Ge\\R-rA,d\Dx),  (3) 

where  the  depth  of  the  distribution  is  |  R-rs  |=  d  (1=0  mm  in  Fig.  1(b)).  To  illustrate  how 

fluorescence  spectra  are  affected  by  varying  the  depth  of  the  fluorophore  distribution, 
Fig.  2(a)  shows  PplX  spectra  computed  with  Eq.  (3)  for  depths  varying  up  to  d= 20  mm 
for  the  absorption  spectrum  labeled  spectrum  1  in  Fig.  3  and  a  constant  value  of 
reduced  scattering  (ps’=1  mm'1).  As  the  object  gets  deeper  in  the  tissue,  the  optical 
properties  demonstrate  an  increasing  influence  on  shaping  the  spectrum  because  of  the 
increasing  pathlengths  traversed  by  light.  The  information  contained  in  the  distorted 
spectra  may  be  distilled  into  a  single  quantity  by  calculating  the  ratio  of  the  signal  at  two 
emission  wavelengths, 


where  the  intensity  values  at  each  wavelength  are  normalized  with  the  relative  signal 
strength  of  an  undistorted  emission  spectrum  (e.g.,  spectrum  labeled  0  mm  in  Fig.  2a) 
implying  that  r  should  be  equal  to  1  for  d=  0  mm.  Normalization  with  the  reference 
spectrum  allows  direct  comparison  between  ratios  computed  with  Eq.  (4)  and 
experimental  data  acquired  with  a  multi-spectral  instrument.  A  critical  aspect  of  Eq.  (4) 
is  that  r  is  independent  of  the  diffuse  excitation  field.  This  means  that,  in  the  point- 
source  approximation  limit,  the  ratio  is  independent  of  the  manner  with  which  the 
imaging  field  is  illuminated.  Therefore,  assuming  that  the  point  approximation  is  valid 
(see  comments  addressing  this  approximation  in  Section  3),  the  equations  presented 
here  are  valid  for  cases  where  the  surgical  filed  is  illuminated  with  a  broad-beam 
source.  Inserting  Eq.  (2)  into  Eq.  (4),  we  find  that  the  logarithm  of  the  fluorescence  ratio 
for  an  infinite  medium  is  linearly  related  to  the  depth  of  the  fluorophores  with  a  slope 
equal  to  the  difference  in  penetration  depth  between  wavelengths  A,i  and  X2, 


lnT 
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1 
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(5) 


This  simple  linear  relationship  suggests  that  depth  can  be  estimated  from  a  simple 
measured  ratio  provided  the  optical  properties  (absorption,  scattering)  of  the  tissue  at 
the  two  wavelengths  are  known.  Also  of  importance  here  is  the  fact  that  this  expression 
provides  a  theoretical  basis  for  interpreting  the  experimental  and  simulation  results 
presented  in  Refs  [7,  8],  However,  the  underlying  photon  diffusion  assumptions  used  to 
derive  Eq.  (5)  restrict  its  application  to  depths  for  which  the  diffusion  approximation 
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applies,  i.e.,  for  objects  located  more  than  1  to  2  mm  underneath  the  surface.  This 
limitation  is  manifest  in  Eq.  (5)  by  considering  the  d= 0  mm  limit,  which  yields  r£^1 . 

In  deriving  Eq.  (5)  two  assumptions  were  made  that  could  potentially  introduce 
errors  in  depth  estimation,  namely  the  medium  was  assumed  to  be  homogeneous  and 
of  infinite  spatial  extent.  While  quantifying  the  effects  of  tissue  heterogeneity  on  depth 
recovery  is  difficult  and  likely  to  vary  for  different  tissue  types,  experience  gained  in 
diffuse  optical  fluorescence  tomography  has  shown  that  using  ratios  of  light  signals  that 
have  sampled  the  same  tissue  volumes  decreases  the  impact  of  heterogeneities  and 
light  coupling  in  and  out  of  the  tissue.  In  other  contexts,  normalization  schemes  have 
been  used  to  bring  tomography  datasets  closer  to  a  form  where  they  can  be  treated  as 
if  they  had  been  obtained  by  interrogating  a  homogenous  medium  [10].  With  respect  to 
boundary  conditions,  the  treatment  for  the  infinite  case  can  easily  be  extended  to  a 
more  realistic  semi-infinite  geometry  which  more  closely  represents  a  surgical  imaging 
configuration,  as  shown  in  Fig.  1(b).  In  this  case,  the  diffusion  equation  is  solved 
ensuring  that  the  fluence  vanishes  on  a  planar  extrapolated  boundary  located  a 
distance  z A  =  2ADa(Ak3  for  air-tissue  interface)  away  from  the  tissue  surface  [11],  The 
corresponding  Green’s  function  is  then 
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where,  for  1=0  mm,  the  distance  between  the  measurement  point  and  the  imaginary 
negative  source  is  given  by\  R-r'm  \=d  +  2z£ .  The  expression  for  the  ratio  of 
fluorescence  becomes 
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where  the  second  term  in  the  brackets  can  be  interpreted  as  a  form  factor  correction  to 
the  expression  obtained  for  an  infinite  medium.  Compared  with  Eq.  (5),  the  semi-infinite 
expression  would  slightly  complicate  the  depth  retrieval  procedure,  since  the  correction 
factor  depends  explicitly  on  the  depth,  which  would  require  the  use  of  iterative  numerical 
resolution  methods.  However,  it  should  be  noted  that  the  form  factor  becomes 
independent  of  depth  for  values  that  are  significantly  larger  than  dc=4ADx.  For  example 
c/c« 2  mm  for  tissue  with  reduced  scattering  values  of  ps’=2  mm'1.  This  implies  that  for 
large  enough  depths  the  linear  relationship  in  Eq.  (5)  is  retrieved  up  to  a  depth- 
independent  scaling  factor  dependent  on  the  tissue  properties  8X  and  D\ 

2.2  Expected  variations  of  measured  fluorescence  ratio  with  depth 

As  evidenced  by  the  expressions  for  the  slope  in  Eq.  (5)  and  the  form  factor  in  Eq.  (7), 
a  critical  ingredient  in  the  accurate  estimation  of  depth  consists  in  the  availability  of  prior 
knowledge  relating  to  tissue  optical  properties  (i.e.,  absorption  and  scattering 
parameters)  at  wavelengths  for  which  fluorescence  signals  are  acquired.  Moreover, 
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initial  assessment  of  whether  or  not  fluorescence  ratio  measurements  can  be  used  for 


depth  retrieval  depends  on  whether  or  not  measurable  signal  changes  with  depth  can 
be  detected  for  optical  properties  that  are  typically  encountered  for  in  vivo  imaging 
applications.  In  order  to  emphasize  this,  Fig.  4(a)  shows  the  fluorescence  ratio  (Xi=650 
nm  and  X2=670  nm)  as  a  function  of  depth  for  the  three  representative  absorption 
spectra  shown  in  Fig.  3.  It  should  be  noted  that  those  spectra  were  computed  based  on 
hyper-spectral  reflectance  data  acquired  in  vivo  with  a  hand-held  fiber-optics  probe  by 
interrogating  human  brain  tissue  as  part  of  a  pilot  clinical  study  involved  12  patients 
undergoing  brain  tumor  resection  procedures  (See  Ref.  [12]  for  details  relating  to  this 
intra-operative  instrument  as  well  as  light  transport-based  algorithms  used  to  retrieve 
the  absorption  and  reduced  scattering  coefficients).  As  shown  in  Fig.  4(a),  the  curve  for 
which  the  difference  in  inverse  penetration  depth  is  the  smallest  (top  curve)  is 
associated  with  a  T-d  relationship  which  is  approximately  linear.  This  can  be  explained 
by  considering  the  linear  limit  of  Eq.  (5),  i.e.,  by  Taylor  expanding  and  assuming  the 
penetration  depth  difference  is  small.  This  finding  is  the  explanation  behind  the 
observation  made  in  Refs.  [7,  8]  that  the  relation  between  r  and  depth  is  linear. 

Equation  (5)  also  indicates  that  the  choice  of  wavelengths  used  in  the  ratio 
affects  the  sensitivity  of  r  to  depth.  Ideally,  the  chosen  wavelengths  should  have 
dramatically  different  optical  properties  to  emphasize  the  spectral  distortion,  which  will 
be  manifested  as  a  steeper  slope  in  the  T-d  relationship.  An  initial  assessment  of  this 
can  be  achieved  by  considering  the  linear  limit  of  Eq.  (5)  showing  that  the  percentage  of 
T-signal  change  per  mm  of  increased  depth  is  approximately 

— ■ -VaV,4  )xioo% i  for//'s  □  jua.  This  implies  that  depth  computation  is 
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facilitated  when  the  difference  in  optical  properties  is  maximized.  However,  under 
similar  light  illumination  conditions,  tissue  that  is  more  absorbing  will  lead  to 
fluorescence  measurements  with  smaller  signal-to-noise  ratio  making  it  more  difficult  to 
detect  at  larger  depths.  As  an  example  of  the  signal  variations  that  can  be  expected  in 
vivo,  optical  property  values  derived  from  the  spectra  shown  in  Fig.  3  can  be  used  to 
estimate  the  expected  levels  of  signal  variations  with  depth.  For  example,  if  ps’=1  mm'1 
and  pax1  =  0.06  mm'1,  absorption  changes  of  Apa  =  pax1  -  Pa2  =  0.001  mm'1,  Apa  =  0.01 
mm"1  and  Apa  =  0.05  mm'1  lead  to  detected  ratio  changes  of  0.4%,  3.7%  and  25.1%  per 
mm,  respectively.  These  variations  further  demonstrate  that  the  choice  of  wavelengths 
used  can  dramatically  affect  the  sensitivity  of  T  to  depth  as  well  as  the  permissible  level 
of  noise  in  the  signal.  Inspection  of  Fig.  3  shows  that  observed  spectral  variations  for 
brain  tissue  should  provide  an  opportunity  to  generate  significant  variations  allowing 
signals  from  different  depths  to  be  discriminated.  Noteworthy  is  the  large  variability  of 
intra-operative  normal  brain  absorption  coefficient  values,  i.e.,  pia  =  0.12  mm"1  to  pa  = 
0.005  mm'1  across  the  spectral  acquisition  domain  ranging  from  630  nm  to  730  nm. 
Further,  to  illustrate  how  fluorescence  spectra  are  affected  by  varying  the  depth  of  the 
fluorophores,  Fig.  2(a)  shows  normalized  simulated  PplX  spectra  computed  with  Eq.  (3) 
for  various  depths  for  the  absorption  spectrum  labeled  spectrum  1  in  Fig.  3  (ps’  =  1  mm' 
1).  At  increasing  depths,  simulations  show  the  influence  of  tissue  optical  properties  on 
the  spectrum  and  its  subsequent  deformation  due  to  the  different  light  paths  traversed. 

An  important  point  that  the  analytical  work  presented  here  is  emphasizing  is  that 
depth  retrieval  based  on  multi-spectral  fluorescence  data  requires  absorption  and 
scattering  properties  of  the  interrogated  tissue  to  be  known  a  priori.  However,  in  a 


13 


context  where  optical  data  is  acquired  in  vivo  during  surgery,  implementing  a  method 
through  which  the  endogenous  optical  properties  are  determined  in  real  time  might  not 
always  be  feasible.  Therefore,  it  is  important  to  understand  how  sensitive  depth 
estimation  is  to  optical  properties  errors  in  order  to  better  assess  the  required  precision. 
As  an  initial  assessment,  the  algorithm  presented  above  was  used  to  retrieve  depth 
using  incorrect  values  of  optical  properties.  To  that  effect,  Fig.  4(b)  shows  curves 
illustrating  the  impact  of  using  the  wrong  optical  properties  when  using  Eq.  (5)  to 
retrieve  depth.  For  the  two  cases  shown  (top  and  bottom  curves),  the  true  optical 
properties  were  assumed  to  match  those  of  spectrum  1  in  Fig.  3.  Using  slopes  and 
intercepts  derived  from  Eq.  (5)  for  Xi=650  nm  and  ^2=670  nm,  values  of  depth  were 
recovered  using  the  incorrect  properties  from  spectrum  2  (top  curve)  and  spectrum  3 
(bottom  curve).  Clearly,  using  the  wrong  optical  properties  can  lead  to  depth  estimation 
errors  of  several  millimeters,  the  magnitude  of  which  increases  with  depth.  The  strong 
reliance  on  the  accurate  tissue  properties  suggests  that  this  approach  should  be 
implemented  with  supplemental  capabilities  to  simultaneously  determine  the  optical 
properties  of  the  tissue  in  the  imaging  field.  In  principle  this  can  be  accomplished  with  a 
system  similar  to  that  shown  in  Fig.  1(a)  where  multi-spectral  or  hyper-spectral  NIR 
transmission  data  is  acquired  and  fitted  for  absorption  and  scattering  spectra  against 
light  transport  solutions  of  the  form  shown  in  Eq.  (6).  A  simpler  approach  would  consist 
in  using  the  fibers-optics  probe  described  in  Ref.  [12]  to  estimate  the  bulk  values  of  the 
tissue  optical  properties  in  the  surgical  cavity  prior  to  imaging  using  the  wide-field  multi- 
spectral  instrument. 
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2.3  Preliminary  experimental  evaluation 

An  initial  experimental  assessment  of  the  proposed  imaging  method  has  been 
performed  by  designing  an  instrument  (see  schematics  in  Fig.  1(a))  for  wide-field 
acquisition  of  multi-spectral  measurements.  Fluorescence  excitation  is  done  using  the 
expanded  beam  of  a  diode  laser  with  wavelength  Xex  =  635  nm  (Power  Technology  Inc., 
Little  Rock,  Arkansas).  Remitted  fluorescence  is  measured  at  five  wavelengths,  X  = 
(675,  685,  695,  705,  715)  nm,  using  10  nm-wide  band-pass  interference  filters  (Omega 
Optical,  Brattleboro,  Vermont)  inserted  into  a  motorized  wheel  coupled  to  an  EXi  Aqua 
charge-coupled  device  (CCD)  (Qlmaging,  Surrey,  British  Columbia).  The  camera  is 
cooled  at  0°C  and  has  a  full  well-depth  of  18,000  e-,  a  14  bits  digitizer,  a  dark  current  of 
0.15  e-/pixel/seconds  and  a  read-noise  of  6.5  e-  when  operated  at  10MHz  readout 
frequency. 

The  proposed  methodology  was  then  evaluated  by  imaging  a  tissue-simulating 
blood  phantom.  A  container  (dimensions:  150  mm  x  75  mm  x  61  mm)  was  filled  with 
water  containing  a  volume  fraction  of  1%  intralipid  (Baxter  Healthcare  Corp.,  Deerfield, 
Illinois)  and  2%  porcine  blood.  An  optically  transparent  cylindrical  inclusion  (diameter: 
5.5  mm,  height:  10  mm)  was  filled  with  20  pg/ml  of  PplX  (Frontier  Scientific,  Logan, 
Utah)  dissolved  in  dimethyl  sulfoxide  (DMSO).  A  moving  stage  allowed  the  inclusion  to 
be  immersed  at  different  depths.  Fig.  5(a)  shows  the  optical  properties  of  the  medium 
computed  using  the  fiber-optics  probe  described  in  Ref.  [12].  The  medium  was 
illuminated  using  an  expanded  laser  beam  (diameter:  13  mm,  P  =  5.5  mW,  A.ex  =  635 
nm).  Fluorescence  images  at  five  different  wavelengths  were  acquired  (d  =  1  to  15  mm). 
A  PplX  spectrum  was  acquired  to  be  used  as  a  reference  that  has  not  been  corrupted 
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by  light  propagation  in  the  medium.  Associated  with  each  fluorescence  image,  a  pre¬ 
contrast  image  (inclusion  filled  with  DMSO  only)  was  acquired  and  subtracted  from  the 
corresponding  image  to  account  for  excitation  light  filter  bleed-through.  For  all  depths, 
the  characteristic  shape  of  the  NIR  PplX  shoulder  was  retrieved  up  to  deformations 
caused  by  diffusion,  as  shown  in  Fig.  2(b).  Fig.  5(b)  shows  the  measured  ratios  for  two 
sets  of  wavelengths:  (A-i,  X2)  =  (675  nm,  685  nm)  (circles)  and  (Xi,  X2)  =  (675  nm,  695 
nm)  (squares).  The  dotted  lines  correspond  to  simulated  ratios  computed  with  Eq.  (6) 
using  the  optical  properties  experimentally  retrieved  with  the  fiber-optics  probe  (Fig. 
5(a)).  Qualitative  inspection  of  the  data  suggests  a  reasonable  match  between  theory 
and  experiment,  although  a  significant  level  of  noise  is  observed  mainly  due  to 
excitation  light  filter  bleed-through,  limited  light  power  delivered  to  the  tissue  surface 
and  variations  of  the  optical  properties  of  the  blood  phantom  during  the  experiment  (as 
much  as  20%  for  the  absorption  coefficient).  Using  a  (non-parametric)  Spearman’s  rank 
correlation  test  a  significant  correlation  was  found  between  the  experimental  data  and 
the  theoretical  predictions.  The  correlation  coefficients  and  P-values  obtained  were  rs  = 
0.90,  P  =  0.0046  (two-tailed)  for  (Xi,  X2)  =  (675  nm,  685  nm)  and  rs  =  0.83,  P  =  0.015 
(two-tailed)  for  (A-i,  X2)  =  (675  nm,  695  nm).  A  linear  regression  was  performed  on  the 
data  with  goodness-of-fit  parameters  r2  =  0.84  and  r2  =  0.64  for  (X-i,  X2)  =  (675  nm,  685 
nm)  and  (Xi,  X2)  =  (675  nm,  695  nm),  respectively.  Another  statistical  analysis  (F-test) 
was  performed  to  test  for  the  presence  of  a  non-zero  slope  in  the  relationship  between 
r  and  d.  We  found  that  for  both  sets  of  wavelengths  the  data  was  significantly  different 
from  a  horizontal  line  (P<0.05). 
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4  Discussion 


Current  methodology  for  brain  tumor  resection  involves  the  use  of  pre-operative 
magnetic  resonance  (MR)  images  co-registered  with  the  surgical  field,  which  is  viewed 
by  the  surgeon  through  a  conventional  surgical  microscope.  More  recently,  intra¬ 
operative  fluorescence  microscopes  have  been  used  providing  surgeons  with  more 
accurate  tumor  detection  capabilities  when  used  in  conjunction  with  the  precursor 
molecule  ALA.  However,  for  a  large  number  of  cases  it  is  suspected  that  even  when  a 
lesion  has  been  fully  resected  based  on  MR  contrast,  trained  visual  inspection,  and 
visible  fluorescence,  residual  tumor  remnants  are  left  behind  underneath  the  surface  of 
the  cavity.  Any  information  that  is  provided  beyond  what  is  currently  seen  using  state-of- 
the-art  fluorescence  surgical  microscopes  is  therefore  likely  to  be  perceived  as  a  major 
advance  in  neurosurgery.  Given  the  current  state-of-the-art  for  fluorescence-guided 
resection  of  brain  tumors  (using  blue  light  excitation)  can  only  detect  on  the  surface  of 
the  cavity,  routine  detection  and  localization  of  tumors  at  1-2  cm  depth  with  NIR  light 
should  lead  to  significant  patient  benefits  in  terms  of  improved  quality  of  life  and 
survival.  Within  the  limits  of  the  approximations  that  were  used  in  the  derivation,  the 
simple  relationships  Eq.  (5)  and  Eq.  (6)  between  fluorophore  depth  and  fluorescence 
ratios  at  two  wavelengths  suggest  that  a  methodology  could  be  implemented  allowing 
intra-operative  estimation  of  residual  tumor  depth  using  non-contact  fluorescence 
imaging  microscopes.  However,  the  relationship  between  d  and  r  depends  explicitly  on 
the  optical  properties  at  the  measurement  wavelengths,  so  a  challenge  for  this  type  of 
approach  is  that  is  must  be  implemented  alongside  an  independent  tissue 
characterization  technique. 


17 


Another  potential  limitation  of  the  approach  is  that  the  point-like  approximation 
has  been  made  in  order  to  derive  a  simple  r-d  relationship  from  which  depth  could  be 
retrieved  with  negligible  computational  burden  when  compared  to,  e.g.,  the  resolution  of 
an  inverse  problem  in  diffuse  fluorescence  tomography.  However,  the  analysis 
presented  in  section  2.3  supports  the  validity  of  this  approximation  by  showing  that 
there  is  a  strong  correlation  between  the  theoretical  model  and  experimental 
fluorescence  ratio  data.  It  should  be  clear  though  that  depth  retrieval  based  on  this 
approximation  will  not  be  exact  because  the  distribution  of  fluorescence  emanating  from 
tumors  is  clearly  not  point-like.  In  reality,  the  assumption  that  is  made  is  that  most  of  the 
fluorescence  that  is  directly  detected  over  the  surface  of  tissue  (1=0  mm  in  Figure  1(b)) 
will  emanate  from  molecular  sources  directly  below  the  surface  area  where 
fluorescence  is  detected.  As  a  consequence,  on  average  the  detected  signal  should  be 
associated  mostly  with  fluorophores  in  the  tumor  that  are  closest  to  the  detection  point. 
In  reality  of  course,  other  parts  of  the  tumor  will  contribute,  so  the  depth  that  is  actually 
retrieved  from  the  signal  will  provide  an  approximate  value,  which  in  effect  will  be  a 
weighted  average  of  the  depth  from  all  contributing  tumor  areas.  The  required  level  of 
precision  for  fluorescence-guided  surgery  (as  attested  by  contributing  neurosurgeons)  is 
such  that  the  precision  reached  for  depth  prediction  with  fluorescence  ratio 
measurements  can  be  a  few  millimeters  off.  Full  validation  of  the  precision  of  the  point¬ 
like  approximation  for  depth  retrieval  will  be  achieved  in  future  studies  where  an 
optimized  instrument  will  be  used  in  the  scope  of  in  vivo  studies.  It  should  also  be  noted 
that  the  point-like  approximation  has  been  used  in  published  work  before,  such  as  in 
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Ref.  [13]  where  time-domain  data  has  been  used  to  predict  the  depth  of  a  fluorescent 
inclusion  of  finite  spatial  extent. 

Because  NIR  light  is  exponentially  attenuated  in  tissue  (modified  Beer-Lambert 
law),  the  information  conveyed  by  epi-illumination  imaging  is  surface  weighted. 
Therefore,  contributions  coming  from  the  auto-fluorescence  of  tissue  can  significantly 
diminish  the  ability  to  image  deeply  buried  sources  of  fluorescence  as  well  as 
significantly  compromise  the  validity  of  the  point-like  approximation.  However,  the 
majority  of  auto-fluorescent  endogenous  molecules  have  excitation  spectra  in  the  visible 
part  of  the  electromagnetic  spectrum.  Since  the  illumination  with  the  proposed 
methodology  is  in  the  far-red  part  of  the  spectrum  (A,ex  =  635  nm),  auto-fluorescence  is 
expected  to  be  small  when  the  method  is  used  in  vivo.  Moreover,  an  oftentimes 
dominating  source  of  tissue  auto-fluorescence  is  the  skin.  However,  in  the  brain  tumor 
resection  application  that  is  presented  in  this  manuscript,  skin  is  not  an  issue  because 
fluorescence  imaging  is  performed  intra-operatively  after  cranium  opening  [14].  Perhaps 
more  convincing  in  addressing  the  issue  of  auto-fluroescence  however,  is  the  fact  that 
the  level  of  ALA-induced  PplX  fluorescence  from  brain  tumors  (including  high-  and  low- 
grade  gliomas,  meningiomas  and  metastasis)  is  two  to  three  orders  of  magnitude  larger 
when  compared  to  normal  adjacent  brain  tissue.  The  high  specificity  of  this  contrast 
mechanism  is  such  that  fluorescence  emanating  from  areas  other  than  the  tumor  is  in 
the  vast  majority  of  cases  negligible  [1 5], 

Therefore  it  is  expected  that  the  validity  of  the  point-like  approximation  will  not  be 
compromised  by  the  presence  of  auto-fluorescence  and  non-specific  fluorescence  for 
such  clinical  applications  as  the  resection  of  brain  tumors  based  on  ALA-induced  PplX 
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fluorescence.  Moreover,  the  preliminary  results  presented  in  this  work  show  that 
fluorophore  depth  correlates  with  measured  fluorescence  ratios  despite  the  fact  several 
approximations  were  made  in  deriving  the  analytic  light  transport  solutions.  This 
provides  strong  evidence  that  a  non-contact  broad-beam  fluorescence  imaging 
approach  can  be  implemented  and  used  intra-operatively  in  order  to  estimate  the  depth 
of  fluorescence  lesions  in  a  manner  that  is  minimally  disruptive  to  the  surgical  workflow. 
However,  using  a  depth  retrieval  algorithm  based  on  Eq.  (6)  with  the  fluorescence-ratio 
data  presented  in  Fig.  5(b)  leads  to  depths  errors  which  can  only  be  minimized  through 
instrument  modifications  reducing  the  level  of  noise  in  the  data.  Future  work  will  be 
concerned  with  those  issues  to  be  addressed  by  developing  an  imager  using  a  light 
source  with  increased  power  as  well  as  light  filtering  components  minimizing  excitation 
light  filter  bleed-through.  Significant  improvements  in  signal-to-noise  characteristics  for 
the  fluorescence  ratios  are  expected  allowing  depth  to  be  retrieved  from  the 
measurements.  Further  validation  of  the  approach  will  be  conducted  by  integrating 
multi-spectral  measurements  into  the  neurosurgical  tumor  resection  workflow  allowing 
for  an  actual  assessment  of  the  in  vivo  level  of  signal  available  as  well  as  the 
development  of  algorithms  optimizing  the  choice  of  wavelengths  insuring  maximization 
of  fluorescence-ratio  variations  with  depth.  Estimation  of  tissue  optical  properties 
(absorption  and  scattering  at  the  measured  wavelengths)  will  be  achieved  by  using  the 
hand-held  fiber-optics  probe  that  is  currently  used  in  the  scope  of  a  brain  tumor 
resection  clinical  study.  The  work  presented  in  this  letter  is  the  first  step  towards  the 
development  of  a  method  that  can  be  readily  applied  to  localize  sub-surface  pathologies 
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during  wide-field  intra-surgical  fluorescence  imaging,  facilitating  more  complete  tumor 
resections. 
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Fig.  1  (a)  Non-contact  multi-spectral  imaging  instrument  used  for  fluorescence-ratio  detection,  (b)  Light 
transport  is  modeled  assuming  the  cavity  can  be  approximated  as  a  semi-infinite  diffusive  medium  with 
homogeneous  optical  properties. 


(a)  (b) 

Fig.  2  (a)  Simulation  results  showing  PplX  spectra  for  molecules  located  at  four  different  tissue  depths, 
(b)  Experimental  spectra  measured  for  a  tissue-simulating  blood  phantom.  The  PplX  spectra  shown  are 
for  different  depths  ranging  from  d= 5  mm  to  c/=15mm.  The  curve  labeled  0  mm  is  the  reference  spectrum 
acquired  with  a  multi-spectral  non-contact  instrument. 
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Fig.  3  Absorption  spectra  acquired  in  vivo  for  normal  human  brain  tissue  using  a  hyper-spectral 
hand-held  fiber-optics  probe. 


(a)  (b) 


Fig.  4  (a)  r  as  a  function  of  depth  for  the  three  absorption  spectra  shown  in  Fig.  3,  (b)  Expected  versus 
computed  depth  for  different  absorption  spectra.  The  curves  labeled  spectrum  2  and  spectrum  3  were 
computed  using  Eq.  (5)  assuming  that  the  optical  properties  were  those  associated  with  spectrum  1. 
Deviations  with  respect  to  the  darker  curve  ( spectrum  1 )  should  be  interpreted  as  depth  errors. 
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—  705  -  675  nm  (theory) 
• —  685  -  675  nm  (theory) 
•  705-675  nm(meas) 

■  685  -  675  nm(meas) 


(a)  (b) 

Fig.  5  Experimental  data  acquired  using  a  tissue-simulating  blood  phantom:  (e)  Optical  properties 
(absorption  and  reduced  scattering)  measured  with  the  fiber-optics  probe,  (f)  Simulated  (theory) 
versus  measured  (meas)  fluorescence  ratios  for  two  sets  of  wavelengths  showing  correlation  with 
depth. 


25 


